International Journal of Modern Engineering Research (IJMER) 
www.iimer.com Vol. 3, Issue. 6, Nov - Dec. 2013 pp-3537-3543 ISSN: 2249-6645 

Comparison of Conventional and "Satlevel" Collocation Model 
for Production of Contour Map of Yanbu Industrial City 



K. F. Aleem 1 ' 2 H. Alghamdi 1 and M. M. Khayyat 3 

1 Department of Geomatics Engineering Technology, Yanbu Industrial College, Yanbu Industrial City, Saudi Arabia 
Surveying and Geo informatics Programme, Abubakar Tafawa Balewa University, Bauchi Nigeria 
Urban Planning Department, Survey Section, Royal Commission for Jubail and Yanbu 

ABSTRACT: Contour map is very useful as its portrays the three dimensional (3D) configuration of the area under study. 
Data acquisition by conventional method of production of such maps has been very tedious. There is need for simple and 
quick method of production like the use of "Satlevel" Collocation. This is a new method of geoid determination in which the 
ellipsoidal height is combined with orthometric height to model the geoid. Geoid determinations using integration of 
geodetic levelling and Global Navigation Satellite System (GNSS), which enables the geoid to be determined in patches, are 
taking over from the empirical geoid determination. Geoid so determined can be applied with ellipsoidal height to get 
orthometric height which height users always prefer. This can then be used for production of contour map. In this work, 
Conventional method was use to acquire data for orthometric height. Also, " Satlevel" Collocation model was used to 
generate orthometric heights in Yanbu Industrial city. The two data were used to produce the contour map of the study area 
using AutoCAD Civil 3D software. "Satlevel" Collocation model enables production of contour map without going through 
all the processes of data acquisition by levelling. The two results were compared and show that there is no significant 
difference between the new map and the other from the acquired data. This method can be extended to other parts of Saudi 
Arabia. 
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I. INTRODUCTION 

1.1 Contour Map: 

Maps are graphic representations of places that use point, line and area symbols, as well as colour, to show how 
selected human and physical features are located, arranged, distributed, and related to one another on the surface of the 
Earth. The complex nature of our environment and the unlimited demand for the uses of maps have made it difficult to have 
a single map that will satisfy all purposes. Therefore, different types of maps exist. Prominent among them is contour map. A 
contour map is made up of contour lines and depicts the shape of any portion of the earth surface. A contour line is an 
imaginary line drawn on the map to connect points with the same elevation or height on, above or below the earth surface to 
form a continuous line shown in brown colour on the map. Historically, Different types of maps were produced for the 
Kingdom by the foreign agencies, educational institutions and governments. The maps includes but not limited to the 
following: 

Historical city maps of Mecca, Al-Madinah and Riyadh from Western Arabia and the Red Sea were produced by 
the Naval Intelligence Division of the Great Britain in 1946; while those of Dharan, Jeddah and Riyadh were produced by 
the United States of America (U.S) Department of State 1983 and 1986. Furthermore, the shaded relief map of Saudi Arabia 
was produced in 1974, reviewed in 1991 and 2003 by the United States Central Intelligence Agency (CIA). Thematic map 
from CIA Atlas of the Middle East was produced in 1993. Topographical maps which include contour were produced as part 
of Operation Navigational chart series by U.S. Defence Mapping Agency Aerospace Centre. In the recent time, U.S. 
National Imagery and Mapping Agency have produced the topographical maps of Saudi Arabia in a scale of 1:250,000 
[PCL, 2012]. 

Presently, the Ministry of Municipal and Rural Affairs (MOMRA) is currently handling mapping and mapping related 
activities in the Kingdom. [Baris, 2012] Apart from production of maps in different scales, different maps are produced 
depending on the need and purposes by other agencies such as the General Commission for Survey (GCS). King Abdulaziz 
City for Science and Technology (KACST), Royal Commission for Jubail and Yanbu, Saudi Geological Survey, Companies 



www.ijmer.com 



3537 I Page 



International Journal of Modern Engineering Research (IJMER) 
www.iimer.com Vol. 3, Issue. 6, Nov - Dec. 2013 pp-3537-3543 ISSN: 2249-6645 

like the Saudi Aramco, Schlumpberger Consulting Engineers Inc., Educational Institutions such as King Saud University, 
King Abdulaziz University, Yanbu Industrial Colleges and others with Geomatics and GIS related disciplines. 
In Yanbu Industrial City, the Urban Planning Department of the Royal Commission for Jubail and Yanbu, Royal 
Commission at Yanbu is currently handling all issues of Survey and Mapping under the Survey section that comprises of 
CAD and GIS sections. The Department produces different types of maps, which include topographical map that contained 
the contours either as a direct labour or contract to reputable mapping companies. 

The contour maps, like most complex ones, are designed using index contours. Contour lines are shown in brown colour 
with the value of height written on it. The height is typical measured with reference to a known Datum. The datum usually 
adopted is the geoid [Aleem et al, 2011 and Williams, 2012] . 

The geoid is the surface which coincides with that surface to which the oceans would conform over the entire earth, 
if free to adjust to the combined effect of the earth's mass attraction (gravitation) and the centrifugal force of the Earth's 
rotation. Specifically, it is an equipotential surface, meaning that it is a surface on which the gravitational potential energy 
has the same value everywhere with respect to gravity. The geoid surface is irregular, but considerably smoother than earth's 
physical surface. Sea level, if undisturbed by tides, currents and weather, would assume a surface equal to the geoid when 
the observation is carried out for a numbers of years, usually 18.61years. [Aleem, 1996, Olaleye et al 2010; Aleem et al 
2011, Aleem, 2013] 

Determination of the geoid has been one of major challenges of geodesists. Gravity data have been used in the past 
with stokes integration and other approaches. These methods are time consuming, expensive and laborious. GNSS provide 
WGSD84 ellipsoidal heights and when compared with orthometric height, from geodetic levelling, allows for the 
computation of the geoid, or the geoid-ellipsoid separation in the region of the survey [Aleem, 1996, Olaleye et al 2010; 
Aleem et al 2011, Aleem, 2013]. These height differences were used to derive the geoidal models called "satlevel" 
collocation model [Aleem, 2013] 

The geoid model will give geoidal undulation at every point of observation. This can be substituted with the 
ellipsoidal height from GNSS observation to get the orthometric height has given in Equation 1.1 [Aleem, 1996, Olaleye et al 
2010; Olaleye et al 2011 Aleem, 2013]: 

H = h - N (1) 

Where: 

H = Orthometric height 
h = Ellipsoidal height 

N = The Geoid-Ellipsoid separation Geoidal undulation 
From the above Equation 1, orthometric height can be computed, which can be plotted on the map as spot heights. The spot 
heights can be interpolated to produce the contours. 

1.2 Interpolation of Contour: 

Traditionally, Interpolations of contours are often done by rough estimate. The procedure is to first plot all the 
available spot heights on the map and estimate the point for the contour. However, a more accurate but time consuming 
method used to be adopted is to use the formula [Aleem, 201 1]: 




Where: 

PC is distance to the point of contour 

D is the distance between two spot heights 

Z is the difference between the two spot heights 

i is the difference between the contour and one of the Spot heights. 
Distance PC will then be scaled and marked on the map, which are later connected to form the contour. The use of 
information Technology in surveying and geoinformatics has made software usage to be part of geomatics. Contour maps are 
nowadays prepared easily with software. 

1.4 Contour Map Software: 

There are different types of software in the market for production of contour, digital terrain and three dimensional 
surface model of any part of the world. This study used AutoCAD Civil 3D software. 

1.4.1 AutoCAD Civil 3D. 

A surface in Civil 3D is built on the basis of mathematical principles of planar geometry. Each face of a surface is 
based on three points within a circum-circle (a circle that passes through each of the vertices of a polygon) forming a triangle 
and defining a plane. Each of these triangular planes shares an edge with another, and a continuous surface is made. This 
methodology is typically referred to as a Triangulated Irregular Network (TIN). On the basis of Delaunay triangulation, this 
means that for any given (x,y) point, there can be only one unique z value within the surface [Golden Software, 2012] 



1.4.2 3DField 

3DField is a contouring surface plotting and 3D data program that runs under Microsoft Windows NT/XP/ Vista/7. 
3DField converts the data into contour maps and surface plots. It creates a 3D map or a contour chart from the scattered 
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points, numerical arrays or other data set. All aspects of 2D or 3D maps can be customized to produce the exact presentation 
desired [Galouchko, 2012]. 

1.4.3 Other Contour Map Software 

Other Contour Map Software includes: Garmin Basemap, Mapviewer, Contour storyteller, Function Grapher, Auto 
plotter, Visual data, Filter test, AutoDEM Li Contour and several other software are available for production of contour map. 
Any of these Software can also be used for the production of Contour map of the study area. 

1.5 Scope and Limitations of the Study:- The research covers the data acquisition, processing and presentation and 
comparison of the final product all were carried out in Yanbu Industrial City while the following are the limitations of this 
research: 

1. The work is only carried out in Yanbu Industrial City. 

2. Automatic level was used to obtain orthometric height as against the use of digital or geodetic level. 

1.6 The study Area 

The study was carried out in Yanbu Industrial City (popularly known as 'Yanbu Al-Sina'iya' in Arabic which 
literarily means Industrial Yanbu) in Madina Province of Saudi Arabia. The city was established around 1975, located on the 
Coast of Red Sea about 350km North of Jeddah. With Latitude 23°59'57.840"N (23.9994) and longitude 38°13'39.000"E 
(38.2275) 




Figure 1: Map of Saudi Arabia and Location of Yanbu Industrial city 

Yanbu Industrial City has a section called the Royal Commission (RC). Its residents are mostly expatriates mostly 
from Europe, America, and other Western Countries. RC Yanbu Existing Area is about 185 km 2 and the expansion Area is 
about 420 km 2 . The industrial developed land is estimated to be 3240 Hectares [RCJY, 2012]. The presence of industries 
such as the oil refineries, plastics facility and several other petrochemical plants in the existing and expansion of Yanbu 
Industrial City will require detailed map and updating of the existing map of which contours have a prominent role to play. 
Contour maps will also enhance monitoring of the existing structures. Therefore, a quick and fast method such the use of 
"Satlevel" collocation model will minimise the rigour of data collection of the conventional method. 

II. AIMS AND OBJECTIVES 

The aim of this work is to produce contour map of the study area from data generated using conventional levelling and 
"Satlevel" Collocation Model for comparison. 
The objectives are: 

1 . To acquire data for ellipsoidal and orthometric heights of the study area. 

2. To compute geoidal coefficients of the "Satlevel" Collocation Model in the study area using Least squares adjustment. 

3. To generate data for production of contour map. 

4. To plot the contour map of the study area using AutoCAD Civil 3D software for the two methods 

5. To compare the contour maps from the two methods. 

III. METHODOLOGY 

3.1 Materials and Data 

The equipment needed for the exercise are: 

GNSS receiver and its accessories to acquire data for ellipsoidal height 
Digital Level and its accessories to acquire data for orthometric height 
Computer and its accessories for computation and analysis. 
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AutoCAD Civil 3D software for plotting the contour 
Any software or program that can implement least square adjustment. 

In this exercise Microsoft Office excels was used for all computations, "Orthometric Height on Fly" developed by Aleem 
[2013] was used to validate the results. 



3.2 "SatleveF Collocation Model: 

Detailed explanation of the derivation and assumption of "Satlevel" Collocation model can be seen in (Aleem, 2013). The 
model is of the form: 

N i = N L + A x (cos 3 cos >l. 4- sin 2 cos <p t cos A t + cos 3 A t cos + sin 2 A t cos cos + 

cos sin A t + sin <p t cos <p t sin A t + cos A t cos <p t sin A t + sin A t cos ^ J + 

cos <p t sin ^. + sin + cos A t sin ^. + sin A t sin ^. J + r t 

[Aleem, 2013] 

Where: 

N L is the long wavelength part of the geoid undulation in the area. 

Ai, A 2 and A 3 are the geoidal coefficients which are unknown coefficients to be determined. 

□ and X are the WGS '84 geodetic coordinates (Latitudes and Longitudes) 

Y[ is residue at an observation point. 
The "satlevel" collocation model geoidal coefficients of the area were computed using least squares adjustment observation 
equation method. 

3.3 Data Acquisition: Levelling operation was carried out to obtain data for orthometric height. The study area is along the 
coast, therefore orthometric correction was neglected. 

GNSS observation was carried out to acquire data for geodetic coordinates which includes the geodetic latitude, geodetic 
longitude and ellipsoidal height. Other relevant data such constant for the referenced ellipsoid were collected from various 
literature and INTERNET websites for data analysis and processing. 

3.4 Data Processing: 

Levelling reduction was done; the reduced levels were the orthometric height of each of the points. GNSS 
observation was processed to get the three dimensional coordinates. These are the geodetic latitude (□), geodetic longitude 
(X) and ellipsoidal height (h), which were used in "Satevel" collocation model. 

The "Satlevel" collocation model geoidal coefficients were computed and used to obtain the geoidal undulation of each point 
in the area, which were connected to get the geoidal surface. 

The geoidal undulations and Ellipsoidal heights obtained from the result of GPS observations were used in Equation 1 to get 
the orthometric height of each of the points in the study area. 

The orthometric heights generated from "Satlevel" collocation and the geodetic coordinates of the area were plotted into 
contour map using AutoCAD Civil 3D software (Figure 5). 

The orthometric heights acquired using geodetic levelling method and the geodetic coordinates of the area were also plotted 
into contour map using AutoCAD Civil 3D software (Figure 6). 

The fallowing procedure was used in producing contour maps in Auto CAD Civil 3D 2013: 
- The data was imported as text file as CSV (Comma delimited) to create the points of the surface. 

The surface was created by a Triangulated Irregular Network (TIN) method with lm interval. 
-The major contours lines were labelled and brown coloured was applied for the contour to follow the convention. 



IV. RESULTS AND DISCUSSIONS 

4.1 Results 

The some of the data used for this research acquired in Yanbu Industrial City were randomly selected and used to plot chart 
(Figure 2): 
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Figure 2: Orthometric and Ellipsoidal Heights of part of Yanbu Industrial City 
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The values of the geoidal coefficients computed using least square adjustment observation equation as: 
N L =7.85995128; A^O.36447485; A 2 =-0.5309024; A 3 =0.15167311 
The results of the "Satlevel" collocation geoidal undulations were plotted into chart (Figure 3): 
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Figure 3: Geoidal Surface of parts of Yanbu Industrial City 




Figure 4: Orthometric Heights of parts of Yanbu Industrial City 



The result of the "satlevel" collocation orthometric height were plotted into chart (Figure 5): 
The results of the Orthometric from conventional levelling were plotted into chart (Figure 6) 



The "Satlevel" collocation contour map (Figure 6): 
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The contour map from levelling method (Figure 5) 




Figure 6: Contour Map from Levelling of parts of Yanbu Industrial City 
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4.2 Discussion 

In this chart shown in Figure 2, ellipsoidal and orthometric height follows the same pattern is an indication that the 
two are true representation of the same terrain. 

Figure 4 shows the geoidal surface of Yanbu Industrial City. In this chart, the geoidal Undulations of 3 of the points were 
outrageous, which is an indication that there was presence of outliers in the observations. These observations should have 
been removed because they are more than 3 standard deviations [Heliani et al, 2004] . This might have created impact on the 
map shown in Figure 6. When the outliers were removed, the chart produced a continuous and smooth surface. This is an 
indication that "Satlevel" Collocation like any other analyses required removal of outliers from the observation. 

V. CONCLUSION AND RECOMMENDATION 

5.1 Conclusion 

In this study, levelled heights were established along with GPS observation in some parts of Yanbu Industrial city 
to model the geoid in the study area. We have coordinated some of the points collocated with both GNSS and levelling in the 
area. 'Satlevel" Collocation Model was used to get the geoidal coefficients, which were used to get the geoidal undulation for 
each of the points. Ellipsoidal heights from GNSS observation were used to get the orthometric heights, and used to plot the 
contour with AutoCAD Civil 3D software. Orthometric height information were also used to plot another map using the 
same AutoCAD Civil 3D software. The two maps were compared with the each other and shows that there is no significant 
difference between the "satlevel collocation and conventional methods of producing contour map. 

5.2 Recommendations: 

The area of coverage needs to be extended and more data should be acquired to improve the results. 
Contracts for the production of maps should be awarded to the Institutions of Higher learning in the Kingdom. This will 
enable the participating staff to acquire more experience, so as to improve the qualities of graduates for future challenges. 
A local geoid model for Saudi Arabia is hereby recommended to commence. 
The method can be applied in any other place. 
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